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ABSTRACT 

We examine the prospect of using the observed abundance of weak gravitational lenses 
, to constrain the equation-of-state parameter w = p/p of the dark energy. Dark energy 

modifies the distance-redshift relation, the amplitude of the matter power spectrum, 
and the rate of structure growth. As a result, it affects the efficiency with which 
1"**- ! dark-matter concentrations produce detectable weak-lensing signals. Here we solve 

the spherical-collapse model with dark energy, clarifying some ambiguities found in 
the literature. We also provide fitting formulas for the non-linear overdensity at virial- 
ization and the linear-theory overdensity at collapse. We then compute the variation 
in the predicted weak-lens abundance with w. We find that the predicted redshift dis- 
tribution and number count of weak lenses are highly degenerate in w and the present 
matter density f^o- If we fix f^o the number count of weak lenses for w = —2/3 is a 
factor of ~ 2 smaller than for the ACDM model w = — 1. However, if we allow fio 
£NJ , to vary with w such that the amplitude of the matter power spectrum as measured 

' by the Cosmic Background Explorer (COBE) matches that obtained from the X-ray 

cluster abundance, the decrease in the predicted lens abundance is less than 25% for 
— 1 < w < —0.4. We show that a more promising method for constraining the dark 
energy — one that is largely unaffected by the Qq — w degeneracy as well as uncertain- 
ties in observational noise — is to compare the relative abundance of virialized X-ray 
lensing clusters with the abundance of non-virialized, X-ray under luminous, lcnsing 
halos. For aperture sizes of ~ 15 arcmin, the predicted ratio of the non-virialized 
to virialized lenses is greater than 40% and varies by ~ 20% between w = — 1 and 
w = —0.6. Overall, we find that if all other weak lensing parameters are fixed, a survey 
must cover at least ~ 40 square degrees in order for the weak lens number count to 
differentiate a ACDM cosmology from a dark-energy model with w = —0.9 at the 3cr 
■ level. If, on the other hand, we take into account uncertainties in the lensing param- 

eters, then the non-virialized lens fraction provides the most robust constraint on w, 
requiring ~ 50 square degrees of sky coverage in order to differentiate a ACDM model 
from a w = —0.6 model to 3c 

Key words: galaxy clusters — weak gravitational lensing — cosmology 
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1 INTRODUCTION 

Observations of distant type la supernovae (SNIa) indicate 
that the universe is undergoing a phase of accelerated expan- 
sion (Perlmutter et al. 1999, Riess 1998). This, combined 
with the flat geometry favored by the cosmic microwave 
background (CMB) measurements (Miller et al. 1999, de 
Bernardis et al. 2002, Halverson et al. 2002, Sievers et al. 
2002, Lee et al. 2001) and the evidence for a low matter- 
density with tt ~ 0.3 (Peacock 2001, Percival et al. 2001), 
suggests that the bulk of the total energy density of the 
universe is in the form of some exotic dark energy with a 
negative equation of state. One of the primary objectives of 



cosmology today is to uncover the origin and nature of this 
dark energy. 

A possible candidate for the dark energy is a cosmolog- 
ical constant A, with an equation of state w — p/p (where 
p is the pressure and p is the energy density of the dark 
energy) strictly equal to —1. Another possibility, and one 
that may find favor from a particle-physics point of view, 
is a dynamical scalar field, termed quintessence, Q. Unlike 
the cosmological constant, the Q-component is both time- 
dependent and spatially inhomogeneous with an equation of 
state w > — 1 that is likely to be redshift dependent. Deter- 
mining the value of w and how it changes with time are key 
to constraining the nature of the dark energy. 
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While the accelerating expansion implies only that w < 
— 1/3, combinations of CMB data, SNIa data, and large- 
scale-structure data suggest that w is most likely in the 
range — 1 < w < —0.6 (Wang et al. 2000, Huterer & Turner 
2001, Bean & Melchiorri 2002, Baccigalupi et al. 2002). 
Though combining these different data sets have provided 
some constraint on w, how w should vary with redshift is 
largely unknown. Particle physics offer several possible func- 
tional forms for the quintessence field's potential V(Q) and 
hence possible scenarios for the time history of to. Nonethe- 
less, determining w's redshift evolution observationally is 
likely to be very challenging (Barger & Marfatia 2001, Maor 
et al. 2001, Weller & Albrecht 2001). 

Strengthening the measured constraint on w and per- 
haps excluding the cosmological constant as the source of the 
dark energy appear, however, to be attainable goals within 
the near future. Since the dark-energy dynamics influences 
both the evolution of the background cosmology and the 
growth of structure, it directly affects many observables. Its 
modification of the angular-diameter distance, the luminos- 
ity distance, and the amplitude of the matter power spec- 
trum, are the primary sources of dark-energy constraint in 
measurements of CMB anisotropics, SNIa, and local cluster 
abundances, respectively. 

In this paper we consider another possible means of con- 
straining w: measurement of weak gravitational-lens abun- 
dances. Weak lensing — the weak distortion of background- 
galaxy images due to the deep gravitational potential of an 
intervening overdensity — provides a powerful technique for 
mapping the distribution of matter in the universe (see re- 
views by Bartelmann & Schneider 2001, Mellier 1999). Here 
we study the impact of the dark energy on the predicted red- 
shift distribution and sky density of weak lenses. Dark en- 
ergy affects the abundance of weak lenses by not only mod- 
ifying the distance-redshift relation and the matter power 
spectrum but also by altering the rate of structure growth. 
In particular, the larger w is the faster and earlier objects 
collapse. An interesting consequences of this is that if we sep- 
arate weak lenses into the two observational classes — those 
that have collapsed and reached virial equilibrium and are 
therefore X-ray luminous and those that are non-virialized 
and hence X-ray underluminous (Weinberg & Kamionkowski 
2002; hereafter WK02) — the abundance of one class evolves 
slightly differently from the other. Therefore the relative 
fraction of these two types of lenses varies with w. This 
observable is especially promising as compared to measure- 
ments of absolute abundances because it is less sensitive to 
uncertainties in both the cosmological parameters and the 
noise in the lensing map. 

This paper is organized as follows. In Section 2 we 
briefly summarize the weak-lensing signal-to-noise estima- 
tor and discuss how we determine the mass- and redshift- 
dependent minimum overdensity required to produce a de- 
tectable weak-lensing signal. Section 3 is devoted to the 
spherical-collapse model in quintessence cosmologies. We 
provide fitting formulas for the non-linear overdensity at 
virialization and the linear-theory density at collapse and 
describe our approach to normalizing the matter power spec- 
trum. In Section 4 we show the resulting effect the dark en- 
ergy has on the weak-lens abundances and in Section 5 we 
present our conclusions. 

Finally, we note that a similar analysis has recently been 



performed by Bartelmann, Perrotta & Baccigalupi (2002), 
although not for the case of non-virialized lenses. Although 
we agree with their general conclusion that the weak lens 
abundance is a potentially sensitive probe of the dark energy, 
our results differ from their results in important details. We 
discuss these differences in Section 4.2. 



2 MINIMUM OVERDENSITY NEEDED TO 
PRODUCE DETECTABLE LENSING 
SIGNAL 

In order to compute the abundance of weak gravitational 
lenses for dark-energy cosmologies we must first determine 
the necessary conditions for a halo of a given density profile 
and redshift to produce a detectable weak-lensing signal. Of 
course the more overdense a halo is relative to the back- 
ground density the more it coherently distorts the nearby 
background galaxies and hence the stronger its lensing sig- 
nal. The detectability of this signal is hampered, however, 
by noise in the weak-lensing map, primary of which is the 
intrinsic ellipticity distribution of the background galaxies. 
The goal is therefore to determine the minimum overdensity 
a halo must have such that it produces a sufficiently large 
signal relative to the noise so as to be detectable. A con- 
venient method for computing this minimum overdensity is 
provided by Schneider's (1996) aperture-mass technique. 

Consider a lens at redshift Zd of surface mass density 
£($) within an angular radius For a source at redshift z s 
the convergence k is given by, 

E(tf) „ c 2 D s 



(1) 



S crit ' 4ttG D d D ds ' 

where D d , D B and Da s are the angular-diameter distances 
between the lens and the observer, the source galaxy and the 
observer, and the lens and the source, respectively. Following 
Schneider (1996), define a spatially-filtered mass inside a 
circular aperture of angular radius 0, 



M ap (0) = / d 2 #K(#)U(\#\), 



(2) 



where U (•&) is a continuous weight function that vanishes for 
# > 6. If U({>) is a compensated filter function, 



d00U(0) = 0, 



(3) 



then M a p can be expressed in terms of the tangential com- 
ponent of the observable shear, 7 t , 

M ap (0) = J d 2 tf 7t (tf)Q(|tf|), (4) 
where the function Q is related to U by 

QW = &f dtftfUtf) -U(fi). (5) 

In this paper we use the I = 1 radial filter function from the 
family given in Schneider et al. (1998): 

Um = ^-(l~x 2 )(±-x 2 ), Q^) = ±^ x \l- X % (6) 

where x = $/9. Taking the expectation value over galaxy 
positions and taking into account the redshift distribution 
of source galaxies then gives, 
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M ap (0) 



, <7t)WQ(l^l), 



(7) 



where (7t}($) is the mean tangential shear on a circle of 
angular radius The function (Z), given by, 



dzs pz(z s )Z(z s ; z d ), 



(8) 



where p z (z s ) is the redshift distribution of source galaxies 
and (Scitz & Schneider 1997) 

lini 2s ^oo S cr it {Zd\ Zs) S cr it c 



Z(z s ;z d ) 



(Zd) 



(9) 



allows a source with a known redshift distribution to be 
collapsed onto a single redshift z s satisfying Z{z s ) — (Z) 
(Seitz & Schneider 1997; Bartelmann & Schneider 2001). 
The source-redshift distribution is taken to be, 



Pz{Zs 



exp [-(zs/zo) 11 



(10) 



r(3/p)z$ 

with (5 — 1.5 and mean redshift (z a ) w 1.5zo = 1-2 (cf., 
Smail et al. 1995; Brainerd et al. 1996; Cohen et al. 2000). 
Finally, assuming the ellipticities of different images are un- 
corrected it can be shown (cf., Kruse & Schneider 1999) 
that the dispersion a M (0) of M ap is 



2 7ro- e 

vm(0) = 

n 



dDDQ 2 (i)), 



(11) 



where n is the number density of galaxy images and a e is the 
dispersion in the galaxies' intrinsic ellipticity. In this paper 
we assume n = 30 arcmin~ 2 and a e = 0.2. The signal-to- 
noise ratio S within an aperture radius 8 is then given by, 

M aD 2{Z)^!°ddd{ lt )(d)Q{d) 



0"M 



(12) 



The tangential shear at (7t)(#), depends on the 
amplitude and shape of the lensing halo's density profile. 
Bartelmann (1995) showed that (7*>(i?) = /t(tf) - («)(#), 
where («) ($) is the dimensionless mean surface mass density 
on a circle of radius $ and /t(#) is the dimensionless mean 
surface mass density within a circle of radius In this paper 
we describe the mass density of lensing halos with the uni- 
versal density profile introduced by Navarro, Frenk & White 
(1996; 1997; hereafter NFW). Thus, for an NFW halo at a 
given redshift with a given mass and mean overdensity rel- 
ative to the background (A = (p pe rt)/pb), we can solve for 
the parameters of the profile (i.e., the scale radius and the 
scale density) and obtain an estimate of (7t)(i?). Details of 
how we solve for the NFW-profile parameters are given in 
the Appendix of WK02. With the density profile known we 
can determine, using equation (12), the expected value of S. 
The minimum mean overdensity, A m i n , needed to produce a 
detectable lens is then given by that overdensity for which 
S > Smin- In this paper we assume S m i n = 5 and 8 = 5', 
unless stated otherwise. 



3 SPHERICAL COLLAPSE IN DARK ENERGY 
COSMOLOGIES 

According to the spherical model of gravitational collapse a 
density perturbation with a nonlinear overdensity A corre- 



1000 



100 



10 



1 1 

_ — w = 


i i i i 
-1 _ 


: --- w = 


-2/3 "--■^^ ' 


W - 


-1/3 \f^^ 


i 





1.2 



1.4 



1.6 
5 



1.8 



Figure 1. The nonlinear overdensity as a function of the linear- 
theory overdensity for three constant-ui models. The full solu- 
tion of the spherical-collapse model predicts collapse to an infi- 
nite overdensity as <5 — > 8 C ~ 1.69. According to the smoothing 
scheme, however, once a mass concentration reaches the virializa- 
tion overdensity A v j r (z), its radius remains constant so that the 
overdensity increases in proportion to the decrease in the back- 
ground density. Shown are the smoothing-scheme solutions for 
mass concentrations that reach the virialization overdensity at 
z = 0. 



sponds to a particular position along the linear-theory evo- 
lutionary cycle. Thus the minimum nonlinear overdensity 
A m i n described above corresponds to a minimum linear- 
theory overdensity <5 m i n ; if an object of mass M at redshift z 
has a linear-theory overdensity 5 > 5 m i n = <5 m i n (M, z), then 
it is sufficiently overdense to produce a detectable weak- 
lensing signal. By determining (5 m i n from the computed A m i n 
we can apply the Press-Schechter (1974) theory to calculate 
the number of halos per unit mass and redshift with 5 > 5 m in 
and hence S > S m i n . We can then find the redshift distribu- 
tion and sky density of weak lenses and how these observ- 
ables vary with w. We will show that for a broad range of 
dark-energy cosmologies a substantial fraction of detectable 
weak gravitational lenses have 5 m i n < S c w 1.69, where 8 C 
is the critical density threshold for collapse. Those objects 
with S < <5 C are commonly thought to be density perturba- 
tions that have not yet reached virialization and are there- 
fore expected to have observational properties that are very 
different from typical virialized lensing clusters. 

In this Section, we present the approach used to map the 
minimum nonlinear overdensity A m i n to a minimum linear- 
theory overdensity 5 m i n for quintessence models (QCDM). 
We describe the dynamical equations of gravitational col- 
lapse in QCDM and give fitting formulas for the nonlinear 
overdensity at virialization, A v i T (z), and the critical density 
Sc- We then discuss how we calculate the abundances of weak 
gravitational lenses, both those with 5 < 8 C and those with 
5 > S c . Below we assume a flat cosmology with a Hubble pa- 
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Figure 2. The radial evolution of a density perturbation that is 
collapsing today according to the spherical-collapse model. The 
ordinate gives the radius, r, in units of the turnaround radius, 
rta, and the abscissa gives the time, t, in units of the overdensity 
collapse-time for the ACDM model. As w increases perturbations 
reach turnaround and collapse earlier, although growth is sup- 
pressed earlier as well. The collapse to a singularity predicted 
by the solution of the spherical-collapse model is avoided by the 
smoothing scheme (thick curves) which yields a constant radius 
once the virialized overdensity is reached. 



rameter h — 0.65, a spectral index n s = 1, a baryon density 
Q h h 2 = 0.02, and Q,q = 0.3, unless stated otherwise. 



3.1 Dynamics 

In quintessence the dark energy is a dynamical, time- 
dependent component, Q, with an equation of state 
parametrized by w = Pq/pq, the pressure divided by the 
energy density. The evolution of the energy density with the 
cosmological scale factor goes as pq oc a~ 3 ' 1+IU ', so that for 
w = — 1 the standard cosmological-constant model, ACDM, 
is recovered. Current observational evidence cannot yet rule 
out a w in the range — 1 < w < — 0.5. 

In order to relate a nonlinear overdensity to a linear- 
theory overdensity in QCDM we must first solve for the 
evolution of the overdensity's radius, R, with time. For a 
spherical overdensity patch with uniform matter density 
ppert = 3M/AttR 3 the evolution is described by the momen- 
tum component of the Einstein equations (Wang & Stein- 
hardt 1998; hereafter WS98): 



^ = jj- l/W + (1 + 3w)p Q \ . 



(13) 



As WS98 pointed out, for w ^= — 1 the space curvature fc per t 
inside the overdensity patch is time dependent. Physically, 
this is because the evolution of the energy density in the 
Q component is evolving independently of the change in ra- 
dius of the overdensity patch. As a result, one cannot assume 



Figure 3. The non-linear overdensity at virialization as a func- 
tion of rcdshift for three constant- w models. As w increases A v i r 
increases because overdensities collapse earlier, when the mean 
gas temperature was higher. For all models A v ; r asymptotes to 
the Einstcin-dc Sitter value of 178 at high redshift. 



that within the collapsing overdensity the rate of change of 
the internal energy in the Q-component, uq, equals the rate 
of work done by the Q-component. That is, because dpQ/dt 
is nonzero unless the Q-component is the cosmological con- 
stant, 



duQ d 

~dT = dt^ v) 



+ -PQ 



dV 

dt ' 



(14) 



where V oc R 3 is the volume of the overdensity patch. There- 
fore equation (13) cannot be cast in the form of a first order 
differential equation as is often done when going from an 
acceleration equation to a Friedman-like energy equation. 
Assuming a constant fc pe rt, as was done in the version 1 
preprint of Lokas & Hoffman (2001), yields significantly dif- 
ferent solutions for the evolution of the radius, R(t), and 
hence for A v j r (z) and So- 
li we combine equation (13) with the Friedman equation 
for the background, 

(a\ 2 8ttG . . ,„„. 

(-) =— (P. + PQ), d5) 

and impose the boundary conditions dR/da\ a=a ^ = and 
-R| o=0 = 0, where a ta is the scale factor at turn-around, then 
for a spherical density perturbation with a given A and red- 
shift z, the unique temporal evolution of the overdensity, 
from linearity to nonlinearity, can be solved (cf., Appendix 
A in WS98). We then have a one-to-one map from A(z) to 
8(z), as shown in Figure 1 for the cases w — — 1, —2/3 and 
— 1/3. The map has a mild w dependence, with a given S 
corresponding to a slightly larger A as w increases. This 
is a consequence of the earlier formation of structure in 
QCDM models relative to ACDM models; overdensities col- 
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lapse faster and are therefore more concentrated for w > — 1. 
This point is well illustrated in Figure 2 where we show the 
growth of a spherical perturbation for the same quintessence 
models. As expected, the larger w is, the earlier structures 
reach turnaround and collapse. 

It can be shown that in the limit 8 — > S c the spherical- 
collapse model predicts that the radius, R, of the overdensity 
goes to zero and hence A — > oo. Of course well before reach- 
ing the singular solution an actual overdensity will virialize, 
thereby halting its collapse. To account for this fact we in- 
voke a simple smoothing scheme in which the radius of the 
matter perturbation is constant with time upon reaching the 
virialized overdensity (see Figure 2). We refer the reader to 
WK02 for details of the smoothing method. 

As described in WS98, the value of A v i r (z) for 
quintessence models, needed here in order to implement the 
smoothing scheme, can be obtained via the virial theorem, 
energy conservation, and solving equations (13) and (15) for 
the overdensity at turnaround. In Figure 3 we show the re- 
sulting numerical solution to A v j r (z). We find that an accu- 
rate fitting function to A v i r (z) for —1 < w < —0.3, modeled 
after the approximation given in Kitayama & Suto (1996) 
for a ACDM cosmology, is 

A vir (z) = 18tt 2 [l + ae b (z)] , (16) 
where 



1.7 



= 0.399 - 1.309(H - 1), 
= 0.941 - 0.205(H ' 938 - 1), 



(17) 



and &(z) = l/OnO?) - 1 = (1/fio - 1)(1 + z) Zw . Since 
structures start to form earlier the larger w is, the mean 
gas temperature in collapsing objects is higher in larger- 
w models. As a result, a greater overdensity is required 
in order for such objects to become bound and virialized, 
explaining why A v i r rises with increasing w. Note, how- 
ever, that for A(z) < A v i r (z) the map from nonlinear to 
linear overdensity has a weak dependence on not only w 
but on fio and redshift as well. The critical threshold for 
collapse today 8 c (z = 0) = S c (z) D(0, fio, w)/D(z, fio, w), 
where D(z,fio,tt>) is the linear growth factor (see WS98), 
also has a weak dependence on fio and w, as shown in Fig- 
ure 4. For 0.1 < fio < 1 and — 1 < w < —0.3, we find that an 
accurate fitting function to 8 c (z), also modeled after the ap- 
proximation given in Kitayama & Suto (1996) for a ACDM 
cosmology, is 

S c (z) = 3(12 2 p )2/3 [l + alog 10 n m ( Z )], 

a = 0.353w 4 + 1.044w 3 + 1.128«; 2 

+ 0.555™ + 0.131. (18) 

Incorrectly assuming that fc pC rt is constant, however, yields a 
(5 C (z = 0) with a much stronger dependence on these param- 
eters, with inferred values for fio = 0.3 of S c (z — 0) ~ 1.5 
and ~ 1.0 for w — —2/3 and w — —1/3, respectively (Lokas 
& Hoffman 2001). 



3.2 Abundances 

Since we are interested in computing the abundances of both 
virialized weak lenses and non- virialized weak lenses we con- 
sider two ranges of overdensity in our lens-abundance calcu- 
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Figure 4. The linear-theory critical threshold for collapse, <5 C , at 
z = as a function of fio for three constant-ui models. 5 C does 
not vary significantly over a wide range in w or Qq. 



lations: (1) <5 m i n < 8 < 8 C , the non-virialized lenses, and, (2) 
8 > S c > Smi n , the virialized lenses. As we showed in WK02, 
both lens types are cluster-mass overdensities. However, 
while the virialized lenses are typically virialized clusters 
that form at rare (e.g., > 3a) high-density peaks of a Gaus- 
sian primordial distribution, the non-virialized lenses corre- 
spond to proto-clusters (e.g., 2a — 3a peaks) — mass over- 
densities that have not yet undergone gravitational collapse 
and virialized, but which have begun to break away from the 
cosmological expansion. These proto-clusters should contain 
galaxies and perhaps a few groups that later merge to form 
the cluster (cf., White et al. 2002). The timescale for col- 
lapse of cluster-mass objects is large, and the overdensities 
can be very large even before they have virialized. It should 
therefore not be too surprising that proto-clusters produce a 
weak-lensing signal that resembles that from virialized clus- 
ters. 

Though the lensing signals may be similar, the two lens 
types are expected to have different observational features. 
In particular, since the X-ray luminosity is a very rapidly 
varying function of the virialized mass, the summed X-ray 
emission from a non-virialized lens should be much smaller 
than that from a fully virialized lensing cluster of the same 
mass. In referring to these proto-clusters as "dark", we thus 
mean that they should be X-ray underluminous. Although 
the mass-to-light ratio of these clusters should be compa- 
rable to those for ordinary clusters, since (1) high-redshift 
clusters may be difficult to pick out in galaxy surveys and, 
(2) proto-clusters will typically have a sky density a few 
times smaller than ordinary clusters, it would also not be 
surprising if these dark lenses had no readily apparent corre- 
sponding galaxy overdensity. Observational evidence of such 
dark lenses has been reported in detections by Erben et al. 
(2000), Umetsu & Futamase (2000), Miralles et al. (2002), 
Dahle et al. (2002), and Koopmans et al. (2000), the lat- 
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ter involving a detection through strong, rather than weak, 
lensing. A more detailed discussion of the features that may 
distinguish dark and virialized weak lenses is given in WK02. 

In order to compute the abundances of virialized and 
dark lenses we need to know the probability that an ob- 
ject of a given mass at a given redshift is in one of the 
above mentioned ranges in overdensity. If we assume Gaus- 
sian statistics for the initial linear-theory density field, then 
the probability that an object's overdensity is in the range 
Si < 8 < 82 is 



P(8 1 <S<8 2 )=er{^j 



(19) 



where 'erf is the error function, v — 8 /a, and a — <j(M, z) 
is the rms density fluctuation of an object of mass M at 
redshift z. From Press-Schechter theory, we know that the 
comoving number density of virialized objects (those with 
8 > <5c) of mass M — 47ri? 3 po/3 in the interval dM that 
are at redshift 2 in a Universe with comoving background 
density p is, 



dn 
dM 



(M,z) 



2 pp 8 c (z) 
7T M 2 a(M,z) 



d\na(M, z) 



x cxp 



dlnM 

8Az) 2 



2<r 2 (M, 



(20) 



We can therefore compute the abundance of objects in the 
overdensity range 81 < 8 < 82 by convolving the above mass 
function of virialized objects with P(Si < 8 < S2)/P(S > 
(5 C ). Specifically, the fraction of objects that can lens relative 
to those that are virialized is, for dark lenses, 



/dark (M, 2) = < 



P(8min <8 <S C ) 
P(8 > 8 C ) 

0, 



< S c 



otherwise, 



(21) 



and for virialized lenses, 



f viT {M,z) = { 



( p(s > a min ) 

P(8 > S c ) 



$min ^ <5ci 

otherwise. 



(22) 



As noted in WK02, the lower the mass of the object the 
larger the minimum overdensity needed to produce a de- 
tectable weak-lensing signal. For low enough masses the min- 
imum overdensity becomes so large that both /dark and / v ir 
approach zero, thereby imposing an effective weak-lensing 
mass threshold. Given / and equation (20) we can compute 
the total comoving number density of weak lenses of a par- 
ticular type. Multiplying by the comoving volume element 
dV c / dz d£l(w) then gives the differential number count of 
lensing objects per steradian, per unit redshift interval: 



dN 



dV c 



dz dD. dz dD. 



f(M)^(M)dM. 



(23) 



By integrating over redshift we can then compute the num- 
ber of dark and virialized lenses we expect to see per unit 
area of sky for a given QCDM model. 

3.3 Normalizing the power spectrum 

In equation (23) the volume term and the two terms within 
the integrand are all functions of w. While the predicted 
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Figure 5. The dependence of as on w as obtained using three 
different approaches: fixing Qq = 0.3 and normalizing to the ob- 
served X-ray cluster abundance (dashed line), fixing f2o = 0.3 and 
normalizing to COBE (dotted line), and allowing Qo to vary with 
w such that the cluster abundance constraint matches the COBE 
constraint (solid line). 





a 




w 



Figure 6. The region in the f2o — w plane where the X-ray clus- 
ter abundance constraint of as, at 95% confidence, overlaps the 
COBE constraint of erg. The gray scale gives the corresponding <r s 
values and the solid line shows whercs the central values match. 
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Figure 7. The comoving volume clement as a function of redshift (left panel) and the comoving number density of virializcd objects 
as a function of mass (right panel) for three constant-M) models. Results arc shown for both the COBE normalization of as with fixed 
Qq = 0.3 (thin lines) and the joint cluster abundance-COBE normalization of erg with Qq = ^o( w ) (thick lines). 



abundance of weak lenses will therefore vary with to, the 
degree to which it will vary depends on the shape and nor- 
malization of the power spectrum of density fluctuations. 
In particular, to compute the abundance of weak-lenses we 
need to know a(M,z). 

For the shape of the power spectrum we use the fitting 
formulas given in Ma et al. (1999) for QCDM models with 
the transfer function and shape parameter for ACDM mod- 
els given by Bardeen et al. (1986) and Hu & Sugiyama (1996, 
eqs. [D-28] and [E-12]), respectively. Since the Q-component 
does not cluster on scales less than ~ 100 Mpc (Caldwell, 
Dave & Steinhardt 1998), at the weak-lensing scales the 
shape of the spectrum does not differ significantly from the 
well studied ACDM shape. 

The normalization of the power spectrum, often ex- 
pressed in terms of as, the rms fluctuation today at a scale 
of 8 h~ x Mpc, is not as well constrained as its shape and 
will in general be a function of w. There are two different 
methods commonly used to obtain the normalization: to fix 
it by the observed X-ray cluster abundance or to fix it by the 
CMB large-scale anisotropics observed by the COBE satel- 
lite. Both approaches have comparable uncertainties; the 
cluster abundance constraint on as has a 20% uncertainty 
at the 2a level (WS98) while the COBE constraint has a 7% 
uncertainty at the la level (Bunn & White 1997). To obtain 
an estimate of how as varies with w so that we may, in turn, 
determine how dN/dzdQ varies with w for dark and virial- 
ized lenses, we will consider three possible approaches. The 
first two involve fixing the cosmological parameters (e.g., 
fio, h, fi b , n s ) and using either the cluster-abundance con- 
strained a s (w) or the COBE constrained as(w). For the for- 
mer we will use the fit given in WS98 and for the latter 
the fit given by Ma et al. (1999); see Figure 5. The third 
approach is to allow the cosmological parameters to be free 
parameters and then jointly match the cluster-abundance 



constraint with the COBE constraint so that each gives the 
same as(w). Since measurements of as are most degener- 
ate with fio, we will let fio be the parameter that varies. In 
Figure 6 we show the region in the fio-w plane where the 
X-ray cluster-abundance constraint, at the 95% confidence 
level, overlaps the COBE constraint. The solid curve shows 
where the central values match, with the resulting range in 
fio (0.3 < fi < 0.4 for -1 < w < -0.4) within observational 
uncertainties (Wang et al. 2000). The corresponding as(w) 
curve is shown in Figure 5. As we will show, the predicted 
weak-lens abundances and how they vary with w strongly 
depend on which as(w) normalization approach is chosen. 



4 RESULTS 

We are interested in determining whether the number count 
and redshift distribution of both dark and virialized weak 
lenses have the potential to constrain w. Another possibly 
useful observable for this purpose is the number count of 
dark lenses relative to virialized lenses. Since dark lenses 
are at an earlier stage of their dynamical evolution as com- 
pared to virialized lenses, those cosmologies that favor a 
faster growth of structure (i.e., QCDM models with larger 
w) will, for a given as, have fewer dark lenses and more viri- 
alized lenses. The ratio of the two is therefore expected to 
vary with w. A priori, this latter observable seems particu- 
larly promising. As discussed in WK02, the ratio of dark to 
virialized lenses is not very sensitive to observational noise in 
the weak lensing maps since observational noise equally af- 
fects the detectability of both types of lenses. Contrastingly, 
uncertainties in observational noise will make it difficult to 
constrain w by simply comparing predicted weak-lens num- 
ber counts with observed weak-lens number counts. 

Before presenting how the above observables are mod- 
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Figure 8. The fraction x of objects with overdcnsities in the range Ata < A < A v i r (upper panels) and < A < Ata (lower panels) 
relative to those objects that are virializcd, A > A v ; r , as a function of mass. The left panels correspond to the COBE normalization of <rg 
with S7o = 0.3 and the right panels correspond to the X-ray cluster abundance normalization of <rg with Qrj = 0.3. For each constant-u) 
model we show the fraction \ at z = (bottom curve), 2 = 1/2 (middle curve), and 2 = 1 (top curve). 



ified by the dark energy we first discuss how each of the 
factors that determine the observed abundance are affected 
by changes in w. Doing so provides both physical insight 
into the results and illustrates the calculational procedure 
discussed in the previous Sections. 

4.1 Preliminaries 

As noted above, the predicted abundance of weak lenses will 
vary with w on account of three factors: the comoving vol- 
ume element, the Press-Schechter comoving number density 
of virialized objects, and the value of /dark/vir [equations (21) 
and (22)]. The degree to which each varies depends on the 
chosen ag(w) normalization. As Figure 7 shows, dV c /dzdQ 
decreases monotonically with increasing w for both fixed firj 
and ilo = flo (w) as given by jointly normalizing as to COBE 
and the cluster abundance. However, because the joint nor- 
malization yields a larger fio with w and a less significant 
decline in as for w > — 1 as compared to the COBE normal- 
ization with fio fixed, the former approach predicts a nearly 
constant virialized object number density with increasing w 
while the latter predicts a significant decrease in the number 
density. 

A similar trend is seen in the functions / v ir and /dark, 



as Figure 8 demonstrates. Here we plot the fraction of ob- 
jects that have not yet reached turnaround (0 < A < A ta ) 
and the fraction of objects that are between turnaround and 
virialization (A ta < A < A v i r ) relative to those objects 
that are virialized (A > A v ir). The figure illustrates sev- 
eral key elements of structure formation according to the 
spherical-collapse model for dark-energy cosmologies. First, 
the fraction, x, of objects in both of these lower-overdensity 
ranges increases with mass in accordance with the hierar- 
chical growth of structure. The fraction also increases with 
redshift since objects are collapsing and evolving toward 
virialization. It is also interesting to note that objects with 
A ta < A < Avir evolve more rapidly as compared to ob- 
jects with < A < A ta - This is demonstrated by the fact 
that at z = 1 the fraction of both types of objects is nearly 
the same though by z = there are more objects that have 
not reached turnaround. Furthermore, the larger w is, the 
greater the difference between the rates of evolution. These 
effects are a consequence of the suppression of structure 
growth in cosmologies with dark energy; namely, growth 
slows down earlier for larger w and those objects that are less 
overdense at a given redshift have greater difficulty overcom- 
ing the repulsive effects of the dark energy and collapsing. 
Finally, the plots show how strongly the fraction depends 
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Figure 9. The normalized rcdshift distribution of virialized lenses (thin lines) and non-virialized lenses (thick lines) for three constant- w 
models. The left panel shows results obtained when erg is normalized to COBE with S7o = 0.3 and the right panel when erg is jointly 
normalized to COBE and the X-ray cluster abundance with f2rj = £lo( w )- The peak of the redshift distributions shift toward lower 
redshifts as w increases because <xg decreases with w. The shift in the peaks is less drastic, however, when the joint normalization is 
assumed. 



on the chosen as normalization, with a significant variation 
with w for the COBE normalization and a fairly small varia- 
tion for the cluster-abundance normalization. This, in turn, 
means that the degree to which the functions / v i r and /dark 
vary with w is highly dependent on the assumed normaliza- 
tion approach. 



4.2 Weak lens abundances 

In Figure 9 we show the predicted redshift distribution of 
virialized lenses and dark lenses for three constant w models. 
For the COBE normalized as with fixed Qo the distributions 
show a fairly strong sensitivity to w. As w increases from 
— 1 to —1/3 the peak of the distributions shift toward lower 
redshifts. Although one might expect the trend to be in the 
opposite direction given that structures form faster for larger 
w models, the effect is counteracted by the decrease in as 
with increasing w. That the decrease in as so overwhelms 
any tendency for structure to form faster for w > —1 is 
not surprising given the weak w dependence in the A - 8 
map (Figure 1) and in the function 5 c (z) (Figure 4). Note, 
however, that the shift in the distributions with w becomes 
much less significant if a joint COBE - cluster abundance 
normalization is assumed. Finally, given that dark lenses are 
likely progenitors of virialized clusters, it is not surprising 
that both normalization approaches predict that the dark 
lenses have a larger mean redshift than the virialized lenses. 

To determine how well the weak lens redshift distribu- 
tions can constrain w we generated mock redshift data and 
determined (using the Kolmogorov-Smirnov test) the proba- 



bility of differentiating two different constant-™ models as a 
function of the number of lenses detected. We found that to 
differentiate a ACDM model from both a w — —0.6 model 
and a w — —0.9 model at the 3a level required, on aver- 
age, approximately 200 weak lenses and 2000 weak lenses, 
respectively. As we show below, this corresponds to a survey 
coverage of ~ 15 and ~ 150 square degrees. Note, however, 
that for sufficiently wide surveys systematic uncertainties 
such as mass-redshift selection effects and lens density pro- 
files might dominate the errors. 

By integrating over the redshift distribution we obtain 
the total number of virialized and dark lenses expected per 
square degree on the sky. As Figure 10 shows, the COBE 
normalization with Qo = 0.3 shows a significant decline in 
the number count as w increases. By w = —2/3 the number 
count of both virialized and dark lenses has dropped by a 
factor of two from the ACDM value. The joint normaliza- 
tion, in which we allow Qo to vary with to, predicts a much 
more mild dependence on w with the number count drop- 
ping by only ~ 20% from w = — 1 to w = —2/3 for both 
lens types. Therefore, while the COBE-only normalization 
approach predicts that the sky coverage needed to distin- 
guish the ACDM model from a w = —0.6 model to 3a is 
only ~ 2 degree 2 , the joint approach requires ~ 15 degree 2 . 
Similarly, to distinguish the ACDM model from w — —0.9 
requires ~ 40 degree 2 and ~ 100 degree 2 , respectively. The 
systematic uncertainties affecting absolute sky density mea- 
surements, such as noise in the lensing maps and uncertain- 
ties in the lens density profiles, are expected to add further 
complications. This suggests that it will be very difficult to 
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Figure 10. The total number of virialized lenses (dashed curves) and non-virialized lenses (solid curves) per square degree as a function 
of w. Thin lines correspond to the COBE normalized erg with Qo = 0.3 and thick lines to the joint COBE— cluster abundance normalized 
us with Qo = Q()(w). While the number count drops by a factor of two between ui = — 1 and w = —2/3 for the COBE-only normalization, 
the drop is much less significant for the joint normalization. 



constrain w using just the number count of either virialized 
or dark lenses without, at the very least, a tighter constraint 
on Qo ■ 

We also note that our results do not agree with the re- 
sults found by Bartelmann, Perrotta & Baccigalupi (2002; 
hereafter BPB). They found that from w = — 1 to w m —0.6, 
the number of virialized weak lenses per square degree in- 
creases by nearly a factor of two. The increase is roughly 
linear up to the maximum after which the number count 
declines steeply. In obtaining these results, however, they 
use the formulas for A v j r and <5 C given in Lokas & Hoffman 
(2001) who assume that the space curvature within a collaps- 
ing overdensity patch is time-independent. As we showed in 
Section 3.1, this assumption is invalid for «i/-l and leads 
to incorrect values for A v i r and S c . To confirm that this is the 
source of our differences, we recomputed the number count 
of weak lenses as a function of w using the algorithm de- 
scribed in BPB (which differs somewhat from ours because 
we are interested in separating lenses into virialized and non- 
virialized types). When we assume the incorrect Lokas & 
Hoffman (2001) values for A v ; r and S c we recover the re- 
sults found by BPB; however, when we assume the values 
for A v i r and <5 C predicted by solving the spherical-collapse 
equations of Section 3.1, we obtain results very similar to 
those described in the preceding paragraphs. 



4.3 Fraction of lenses that are dark 

As mentioned above, the number-count ratio of dark to viri- 
alized lenses is an observable that is much less sensitive 
to observational noise than is the redshift distribution and 
number count of weak lenses. Unfortunately, for aperture 
sizes 6 (defined in Section 2) less than 10' in radius the ra- 
tio is fairly constant over a broad range in w, as we show 
in Figure 11. The ratio varies more strongly if the aperture 
size is increased to 15'. In particular, for 9 — 15' there is a 
~ 20% difference between the ACDM model and w = —0.6, 
so that differentiating the two models to a 3u significance 
requires the detection of ~ 600 virialized lenses or equiva- 
lently a sky coverage of ~ 50 degree 2 . Although using the 
non-virialized lens fraction requires large survey coverage for 
modest constraints on w, its principal advantage (in addi- 
tion to being relatively insensitive to observational noise) is 
that it is not very sensitive to the chosen method of normal- 
ization; for any aperture size both the joint normalization 
and the COBE normalization with fixed f2o yield similar 
dependences on w. Therefore, unlike the case for weak-lens 
sky-density or redshift distribution predictions, uncertain- 
ties in erg and f2o do not strongly affect the predicted ratio 
of dark to virialized lenses. Incidentally, although aperture 
sizes greater than ~ 15' yield ratios with even stronger w 
dependences, noise contributions from large-scale structure 
become significant at such large angular distances from the 
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Figure 11. The fraction of lenses that are dark as a function 
of w for aperture sizes 8 = 3, 5, 10 and 15 arcmin. Results are 
shown for both the COBE normalization with f2o = 0.3 (dashed 
lines) and the joint COBE-clustcr abundance normalization with 
Qo = f^o('i') (solid lines). As 8 increases, the fraction of lenses 
that are dark rises significantly. 



Figure 12. The total number of virialized lenses (dashed line) 
and non-virialized lenses (solid line) per square degree as a func- 
tion of the aperture size 8 for a ACDM cosmology. While the 
number count of virialized lenses peaks at 8 = 5 arcmin and 
declines thereafter, the number count of non-virialized lenses in- 
creases almost linearly for 8 > 5 arcmin. 



lens center (Hoesktra 2002). It is therefore not practical to 
make measurements at radii well beyond 15'. 

As an aside, while the ratio of dark to virialized lenses 
does not have a particularly strong w dependence, it does 
have a strong 9 dependence; only ~ 5% of lenses are dark 
when 8 = 3' but ~ 50% when 9 = 15'. In Figure 12 we 
plot the number of virialized and dark lenses as a func- 
tion of 9 for the ACDM model. As 9 increases from 3' to 
15' the sky density of dark lenses increases from zero to 5 
per square degree while the sky density of virialized lenses 
peaks at 9 — 5' and gradually declines for larger aperture 
sizes. Figure 13 explains this trend. For an overdensity of 
mass M = 5 x 10 14 M Q we plot, as a function of redshift, 
# v ir, the projected angular size of the virialization radius, 
and # m ax, the projected angular size of the maximum radius 
that produces a detectable lens (i.e., # m ax — Rmmx(z)/ Dd(z) 
where = 3M/47rA min (z)). For max > 9 vil an over- 

density can be non-virialized and still produce a detectable 
lensing signal (i.e., a dark lens). However, since 9 defines the 
maximum observable angular scale, for sufficiently small 9 
there is no range in redshift such that 9 > # max > 9 v i r , in 
which case non-virialized overdensities cannot produce a de- 
tectable lens. In general we find that the minimum aperture 
size needed to detect dark lenses is ~ 3'. For larger 8, the 
area below max and above 9 v i T has a substantial relative 
increase while the area below # v j r has just a mild relative in- 
crease. After taking into account the fact that the aperture 
mass M ap (8) decreases with increased 9, this translates to 
an increase in the sky density of dark lenses and a decrease 
in the sky density of virialized lenses for 9 > 5'. The fraction 
of lenses that are dark therefore increases with aperture size. 



5 DISCUSSION AND CONCLUSIONS 

We have examined the possibility of using the measured 
abundance of weak gravitational lenses to constrain a prin- 
cipal property of the dark energy, its equation-of-state pa- 
rameter w. Since dark energy modifies both the background 
cosmology of the universe and the growth of structure it 
will necessarily have an effect on the efficiency of weak lens- 
ing. The goal of this paper was to determine the nature and 
strength of the effect. 

The change in the background cosmology with w influ- 
ences the predicted weak lens abundance in essentially three 
ways. First, the size of comoving volume elements shrink 
with increasing to. Second, the distance-redshift relation is 
modified, thereby shifting the location of the lensing- kernel 
maximum (i.e., where the combination of angular diame- 
ter distances DdsDd/ D s peaks). Third, since the evolution 
of the background matter density is modified by the dark 
energy, the density of a given halo relative to the back- 
ground density changes with w. This, in turn, affects the 
strength of a halo's lensing signal; the larger the overdensity 
the stronger the signal. While the volume term is explicitly 
factored into the expression for the weak-lens sky density 
[equation (22)], the latter two effects are incorporated into 
the signal-to-noise estimator for which we use the aperture- 
mass technique introduced by Schneider (1996). 

The change in the growth of structure with w is some- 
what more subtle. The dark energy modifies both the rate 
of structure growth and the amplitude of the matter power 
spectrum. To determine the former we solved the spherical- 
collapse model with dark energy included. Though growth 
occurs more rapidly as w increases, the overall effect on the 
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Figure 13. The observed angular size of 9 v i T (dashed line) and 
Omax (solid line) as a function of redshift for an overdensity of 
mass M = 5 X 10 14 M Q for a ACDM cosmology. If max > $vir 
an overdensity need not be virializcd to produce a detectable 
lensing signal. However, the range in redshift over which a non- 
virialized lens can be detected is limited by the aperture size 8 
(e.g., thin, dotted lines), which defines the maximum observable 
angular scale. For $ < 3 arcmin virtually no dark lenses can be 
detected. 



A - 5 map, needed to relate the minimum overdensity re- 
quired to produce a detectable lens, A m i n , to a correspond- 
ing linear-theory overdensity <5 m in, is fairly small. Similarly, 
the linear-theory overdensity at collapse <5 C does not vary 
much with w. The effect on A v ir is more significant, how- 
ever. As w increases, structures require substantially greater 
overdensities in order to reach virial equilibrium because 
they collapse sooner, when the universe was younger and 
hotter. 

To determine how the power-spectrum amplitude, as, 
varies with w we considered three possible approaches. One 
was to normalize to the X-ray cluster abundance as was 
done in WS98. Another was to normalize to the COBE 
measurements of CMB anisotropics on large angular scales. 
These two approaches predict similar values of as for the 
ACDM model. However, if all cosmological parameters are 
held fixed as w varies, the values of as are no longer in ac- 
cordance. This is because the cluster abundance approach is 
accounting for the earlier forming, and hence hotter, galaxy 
clusters in models with w > — 1. The COBE normaliza- 
tion, on the other hand, is accounting for the increase in 
the Integrated Sachs- Wolfe (ISW) effect as w increases (cf., 
BPB). Given these differing influences, the two approaches 
are not expected to yield the same as when all the cos- 
mological parameters are held fixed to those of the ACDM 
model while w is varied. This suggests a third approach to 
normalizing the power spectrum; namely, let the parameters 
vary with w such that the cluster abundance normalization 
matches the COBE normalization. In practice we accom- 



plished this by letting just Qo vary with w, as it is the pa- 
rameter most degenerate with as . The resulting range in Qo 
for -1 < w< - 0.4 was found to be 0.3 < Q < 0.4 and 
hence within observational uncertainties. Though all three 
normalization approaches predict that as decreases with w, 
the difference in the magnitude of the decrease between the 
approaches is significant. As a result, each predicts substan- 
tially different variations in the weak-lens abundance with 
w. 

Having determined all the dark energy effects, we com- 
puted the redshift distribution and sky density of weak 
lenses as a function of w. As in WK02, we distinguished 
between two classes of lenses, those that have collapsed and 
virialized and those that have not. This distinction is based 
on the expectation that the virialized lenses, being in a re- 
laxed state, are X-ray and/or optically luminous. The non- 
virialized lenses, being as at earlier stage in the overdensity 
evolutionary cycle, are expected to be X-ray underluminous 
because the observed X-ray luminosity function has a steep 
dependence on the total virialized mass within a halo. Fur- 
thermore, though the typical mass of both lens types is 
~ few x 10 14 M Q , the sky density of galaxies within the 
non-virialized lenses is expected to be smaller than in the 
virialized lenses because they have not yet collapsed and 
hence have larger radii (see WK02 for more details). 

We found that the variation in the redshift distribu- 
tion and the sky density of both lens types with w depends 
strongly on the power-spectrum-normalization approach. If 
Qo is fixed and as is normalized to the COBE measurements, 
there is a significant variation in the abundances with w. 
In particular, the sky density of both virialized lenses and 
non-virialized lenses drops by a factor of two from w = — 1 
to w — —2/3. This decline, a result of the significant de- 
crease in as with w, occurs despite the faster formation of 
structure for w > — 1. If, on the other hand, Qo is allowed to 
vary with w such that the COBE normalization matches the 
cluster-abundance normalization, the redshift distributions 
and sky density change very little with w; between w = — 1 
and w = —2/3 the sky density of both lens types varies by 
just ~ 20%. This insubstantial variation is the result of an 
increase in Qo with w and a less significant drop in as with 
w as compared to the COBE normalization with Qo fixed. 
Obtaining a strong constraint on w from the sky density or 
redshift distribution of weak lenses therefore appears to be 
contingent on improved measurements of Qo from indepen- 
dent observations. 

Perhaps more promising is the possibility of utilizing 
the observed ratio of dark lenses to virialized lenses. Unlike 
measurements of the absolute sky density of weak lenses, 
the ratio is not very sensitive to the amount of observational 
noise in the weak-lensing maps since the abundance of both 
dark lenses and virialized lenses are equally affected by noise. 
Similarly, the ratio does not vary significantly over a wide 
range in cosmological parameters so that uncertainties due 
to the Qo — w degeneracy are minimized. We found that 
for aperture sizes of ~ 15' the ratio varies by about 20%, 
dropping from 0.5 to 0.4, between the ACDM model and 
w = —0.6. We also showed that the ratio of dark to virialized 
lenses increases with aperture size, in effect because larger 
apertures enable the detection of the more extended radii of 
the non-virialized lenses. 

Weak lensing has already been shown to be a powerful 
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probe of the matter distribution in the universe (see e.g., 
Bartelmann & Schneider 2001). It also has the potential to 
help constrain the amount and nature of the dark energy. 
Huterer (2002) showed that given reasonable prior infor- 
mation on other cosmological parameters, the weak-lensing 
convergence power spectrum can impose constraints on the 
dark energy comparable to those of upcoming type la su- 
pernova and number-count surveys of galaxies and galaxy 
clusters. Constraining the dark energy from absolute mea- 
surements of weak-lens abundances will likely prove difficult, 
however. The variation in the weak lens sky density with w 
is sufficiently small that modest uncertainties in Qo (and 
observational noise) can mask the effect of the dark energy. 
More auspicious is the possibility of utilizing the relative 
abundance of dark lenses to virialized lenses to constrain w. 
Future weak-lensing projects such as the the VISTA survey, 
the SNAP mission, and LSST (see Tyson et al. 2002 for a 
discussion of its great promise as a probe of dark energy) 
are expected to provide the wide-field surveys needed for 
this technique to be viable. 
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